Phase Transition in the Two-Dimensional Gauge Glass 
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The two-dimensional XY gauge glass, which describes disordered superconducting grains in strong 
magnetic fields, is investigated, with regard to the possibility of a glass transition. We compute the 
glass susceptibility and the correlation function of the system via extensive numerical simulations 
and perform the finite-size scaling analysis. This gives strong evidence for a finite-temperature 
transition, which is expected to be of a novel type. 
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The gauge glass model, which was originally proposed 
as a generalization of the spin-glass model [Q, has at- 
tracted much attention in relation to the vortex-glass 
phase of high-T,, superconductors H. In three dimen- 
sions, the XY gauge glass model is believed to exhibit 
a finite-temperature glass transition [pUj], in agreement 
with experimental evidence for the vortex-glass phase 
at finite temperatures p|. In two dimensions, on the 
other hand, there has been controversy as to the ex- 
istence of a finite-temperature transition. Equilibrium 
studies of several quantities such as the defect-wall en- 
ergy M and the root-mean-square current H] have sug- 
gested the absence of ordering at finite temperatures; this 
appears to be consistent with experiment, where no evi- 
dence for the glassy phase at finite temperatures has been 
found g. Indeed the gauge-glass order parameter has 
been shown analytically to be zero at any finite temper- 
ature [y. However, the helicity modulus computed via 
Monte Carlo simulations indicates a signal of a glassy 
phase at low but finite temperatures |10|] . Dynamical 
simulations also give conflicting results: Whereas earlier 
study of the current-voltage characteristics has shown a 
strong evidence for the possibility of glass ordering at fi- 
nite temperatures |lj] , later one has been interpreted to 
be consistent with the zero-temperature transition [h2|. 
However, the lowest temperature considered in Ref. |l^ ] 
is apparently higher than the estimated transition tem- 
perature Tc ~ 0.15 (in units of the coupling energy J) 
| jll[ . Indeed the recent study of relaxation dynamics has 
indicated a finite transition temperature Tc « 0.22 [ p^ . 
It should be noted here that the absence of an ordered 
phase does not necessarily imply the absence of a phase 
transition, as we will discuss later. We thus believe that 
the presence/absence of a finite-temperature transition in 
two-dimensions as well as its nature is still inconclusive. 
For the resolution of this, careful analysis of the behavior 
should be performed at sufficiently low temperatures. 

In this work we investigate the two-dimensional 
XY gauge glass via extensive numerical calculations, 
with particular attention to the possibility of a finite- 
temperature glass transition. We adopt the equilibration 



test method in Refs. [ [l4|JT5| to obtain the equilibrium 
and to determine the equilibration time, up to the sys- 
tem size L = 48. ^From the obtained equilibrium con- 
figuration, we compute the glass susceptibility and ex- 
amine it by the finite-size scaling analysis [ [l4| . This re- 
veals remarkable divergent behavior at low but non-zero 
temperatures, providing strong evidence for the finite- 
temperature transition. It is discussed how the presence 
of such a finite-temperature transition can be reconciled 
with the results of existing studies. 

We consider the standard XY gauge glass model on a 
square L x L lattice, which is described by the Hamilto- 
nian 
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where J is the coupling energy between nearest- 
neighboring grains, (pi is the phase of the order parameter 
of the grain at site i = {x, y) (a;, j/ = 1, 2, • ■ • , L), and the 
bond angle Aij 's are taken to be quenched random vari- 
ables distributed uniformly on the interval [0,27r). The 
presence of a glass transition in the system can be con- 
veniently described by the divergence of the glass sus- 
ceptibility, given by xg = X)? ^ij "^ith the correlation 
function of the glass order parameter 



G. 
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where [• • •] denotes the disorder average and (• • ■) the 
thermal average. In the limit the distance between the 
two grains i and j becomes large, the correlation function 
reduces to q^ , the square of the Edwards- Anderson glass 
order parameter 



|(e'^.)l 



(3) 



In an infinite system the glass susceptibility is expected 
to display the critical behavior 

XG-iT-Tcy (4) 

with the scaling relation among the exponents 



7 = (2 - r])v, 



(5) 



where rj describes the power law decay of the correlation 
at Tc and v the divergence of the correlation length ^: 
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as the temperature T approaches Tc. In a system of 
size L, the glass susceptibility has the finite-size-scaling 
form liljl5tl 



XG^L'-''x{L'/''iT~Tc)) 



(7) 



with the appropriate scaling function x, which can be ex- 
amined via extensive simulations at various temperatures 
and sizes. 

In the simulation concerning equilibrium properties, 
equilibration is an important issue. We follow Ref. Pq \ 
to have a criterion for equilibration, and consider two 
replicas a and (3 of the system with the same realization 
of disorder. The susceptibility can be calculated from the 
overlap between the two replicas, according to 
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where A^ = L^ is the number of spins and time t is mea- 
sured in units of the Monte Carlo sweep (MCS). Note 
also that the time-dependent four-spin-correlation func- 
tion defined as 



Xcito) = -^ 
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converges to the glass susceptibility xg in the limit 
to — * cxD. Thus the equilibration time r can be estimated 
as the value of to at which the two expressions, Eqs. (J^) 
and (m, give coincident results. The expected behavior 
of the equilibration time r is 



r oc C" oc (T - TcY 
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in an infinite system near Tc, with the dynamical expo- 
nent z; in a finite system at T = Tc, it is expected 



TOCL^ 
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Equations (JlOj) and (|ll| ) naturally lead to the finite-size 
scaling form |16|| 



T^L^f{L^/''{T-Tc)). 



(12) 



To examine the behavior of the glass susceptibility, we 
have performed extensive simulations at several temper- 
atures, ranging from T = 0.2 to T = 1.0. After the 



equilibration time t estimated as above, the glass sus- 
ceptibility has been measured according to Eq. (||) ex- 
cept for that the data have been taken during sufhcicntly 
large time intervals, from three to ten times of the equi- 
libration time. Namely, to has been chosen to be from 
3t to lOr, depending on the size. We have also per- 
formed independent runs with up to 16000 different dis- 
order configurations, over which the disorder average has 
been taken. Since the different disorder realizations give 
statistically independent thermal averages, the statisti- 
cal error can be estimated from the standard deviation 
of the results for different samples. Here the number of 
the disorder realizations to get sufficiently reliable statis- 
tics for the equilibration time turns out much larger than 
that for the glass susceptibility. Thus it has taken much 
more averages to obtain reliable data for the equilibra- 
tion time than those for the glass susceptibility. In this 
way, for example, the equilibration time at T = 0.3 has 
been estimated to vary from (7± 1) x 10'' MCS (L = 4) to 
(2 ± 0.1) X 10^ MCS (L = 24); it increases rapidly as the 
temperature is lowered, reaching at T = 0.25 the value 
(1.5 ±0.4) X 10^ for L = 12. 

Figure [^ shows the behavior of the obtained glass sus- 
ceptibility Xg with the temperature T for system sizes 
L — 4, 6, 8, 12, 24, and 48. In all cases, periodic bound- 
ary conditions have been employed, and the data points 
without error bars indicate that the errors estimated by 
the standard deviations are smaller than the size of the 
symbols. It is indeed observed that, as L is increased, 
the data points apparently approach the dashed line rep- 
resenting (T — Tc)~'^ with Tc and 7 obtained from the 
finite-size-scaling method below. The finite-size scaling 
analysis of the data is displayed in the inset of Fig. H, 
where the data collapse nicely to the finite-size-scaling 
form given by Eq. Im) . The corresponding exponents and 
the transition temperature are obtained: 



?7 = 0.30 ±0.05, 

iy= 1.14 ±0.07, 

rc = 0.22 ±0.02, 



(13) 



which agree with the results of Ref. [|l3[ . The errors have 
been estimated by the size of the region beyond which 
the data do not scale well. The scaling law in Eq. (|^) 
then gives the susceptibility exponent 



7 = 1.93 ±0.08. 



(14) 



We have also tried to fit the data to the finite-size scal- 
ing form with Tc — 0, and observed that the high- 
temperature data (T > 0.35) appear to fit also the Tc — 
scaling, in agreement with Ref. |l2|. At lower temper- 
atures, however, systematic deviation from the scaling 
has been revealed, which becomes conspicuous for large 
sizes (L = 24, 48) |0]. Note that such low-temperature 
regions, which are crucial in differentiating the two seal- 



ings {Tc — and Tc ^ 0), have not been probed prop- 
erly in previous studies concluding a zero-temperature 
transition. Moreover, even for small sizes such zero- 
temperature scaling yields rj very close to zero, which is 
rather unlikely in view of the large ground-state degen- 
eracy expected in the system |l^. We also perform the 
finite-size scaling analysis of the equilibration time. In 
this case, assuming the scaling form in Eq. (p^), we use 
the values of the exponent u and the transition tempera- 
ture Tc in Eq. (|13), and estimate the dynamic exponent 
z. Figured shows that such scaling is indeed reasonable, 
with the dynamic exponent estimated as 



the energy change (in units of J) E 
thus the desired relation 



Ep = 0{L°), and 



2.4 ±0.3. 



(15) 



In contrast, these equilibration time data hardly fit the 
scaling with Tc = 0, displaying marked deviation at 
T<0.4 |l|]. 

We now examine the implications of such a finite- 
temperature transition, and first consider the scaling be- 
havior of the defect- wall energy, which is defined to be the 
difference in the ground-state energy upon changing the 
boundary conditions along one direction from periodic to 
antiperiodic y,g| . The defect- wall energy fluctuates from 
sample to sample with zero mean, and its typical value, 
which may be taken as the average of the absolute value 
over quenched randomness, scales with the system size L 
according to 
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in the asymptotic domain {L — > oo) [ [L8[ . The sign of 
the exponent 9 then determines the presence or absence 
of long-range order at finite temperatures: Whereas for 
6 positive, the system displays rigidity or ordering, the 
negative value of 9, obtained numerically in two dimen- 
sions ||] , implies ubiquity of long- wavelength fluctuations 
and thus the absence of order. Here it should be stressed 
that the absence of (long-range) order does not neces- 
sarily correspond to the absence of a finite-temperature 
transition, as the Berezinskii-Kosterlitz-Thouless (BKT) 
transition and the associated algebraic order present an 
example [^. In fact it is easy to show that in any 
two-dimensional system described by the Hamiltonian 
(|l|), gapless spin- wave excitations prevent 9 from hav- 
ing a positive value: Suppose that the ground-state en- 
ergy in the periodic boundary conditions (PBC), Ep, 
is smaller than that in the antiperiodic boundary con- 
ditions (APBC). ^From the ground-state configuration 
{(j)l } in the PBC, we can make a configuration sat- 
isfying the APBC by rotating the phases according to 
^i ^i + a;7r/L. It is then obvious that the energy E 

of the new configuration is not smaller than the ground- 
state energy Ea in the same APBC. In the limit L -^ oo, 
upon expanding the energy E of the new configuration 
around the ground-state configuration {(j)- }, we obtain 



Ea-Ep<E-Ep = 0{L''), 



(17) 



which shows that 9 cannot be positive. For Ea smaller 
than Ep, similar argument starting from the ground-state 
configuration in the APBC then yields Ep — Ea < 0{L°). 
Accordingly, it may be the case that the negative value of 
9 results from spin- wave excitations, which destroy long- 
range order but maintain criticality. 

To examine this possibility, we have generated ground 
states in the PBC and in the APBC via extensive simu- 
lations for several system sizes (up to L = 24), which has 
confirmed the value of 9 obtained in the existing stud- 
ies on small sizes [g[. We have then carefully compared 
the two corresponding ground-state configurations, one 
in the PBC and the other in the APBC, and found only 
small difference in the total number of vortices. In partic- 
ular the relative difference appears to decrease with the 
system size [|7| , indicating that no (bulk) vortex excita- 
tion is involved. This suggests that the ubiquitous fluctu- 
ations implied by the negative value of 9 are indeed of the 
spin- wave type rather than of the vortex one [ pO[ . Such 
spin-wave excitations should lead to the algebraic decay 
of the correlation function in Eq. (0) and to the vanish- 
ing Edwards- Anderson glass order parameter defined in 
Eq. (0) at all finite temperatures, thus consistent with 
the result of Ref. P|. It is thus strongly suggested that 
the system displays quasi-long-range glass order below 



the finite transition temperature |21|, characterized by 
the algebraic decay of the glass correlation function. 

We have thus computed the glass correlation function 
via large-scale simulations for the system size L = 48 and 
display the obtained behavior in Fig. 0, where the differ- 
ence in the behavior at the two temperatures T = 0.15 
and 0.50 is spectacular: While the dashed line at T = 
0.50 represents the least-square fit to the usual exponen- 
tial decay r~^e~^^^ with 77 = 0.31 and ^ = 1.90, the alge- 
braic decay at T = 0.15 is manifested by the least-square 
fit (represented by the dotted line) to r~^ with 77 = 0.27. 
The decrease of the exponent 77 with temperature, from 
the value 0.30 (at Tc), refiects the criticality of the sys- 
tem below Tc, confirming quasi- long-range glass order at 
low temperatures. The fit of the data at T = 0.15 to the 
exponential-decay form has been tried, only to yield very 
large deviations as well as inconsistent values of param- 
eters 77 and ^. We have also computed the correlation 
length at various temperatures above Tc and found per- 
fect scaling to the form in Eq. (H) with the values of Tc 
and v given by Eq. (O) (see the inset). 

We finally point out that the existence of the finite- 
temperature glass transition is not inconsistent with 
the experimental results. The temperature scale in the 
gauge-glass model may be obtained from the correspon- 
dence between the model temperature T = 1.15 (in units 
of J) and the real temperature T = 12K in experiment 



| jr2[ . Accordingly, the transition temperature estimated 
here corresponds to approximately 2. OK; such a low- 
temperature regime was not probed in Ref. g and it 
would be interesting to investigate the low-temperature 
regime experimentally. 
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FIG. 1. Behavior of the glass susceptibility with tempera- 
ture T. The dashed line is proportional to (T - 0.22)^^'^^. 
Inset shows the fit to the finite-size-scaling formula given by 
Eq. (JtI), with Tc = 0.22, 1/u = 0.88, and 2 - 77 = 1.70. 
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FIG. 2. Fit of the obtained values of the equilibration 
time to the flnite-size-scaling formula given by Eq. (H), with 
Tc = 0.22, 1/u = 0.88, and z = 2.4. 
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FIG. 3. Behavior of the correlation function with distance 
r at two temperatures, T — 0.15 and 0.50. Algebraic decay is 
obvious at T = 0.15. Inset shows correlation length ^ versus 
temperature T — Tc in the log-log scale with Tc = 0.22 ± 0.01. 
The slope of the soild line is fixed to the value —v — —1.14. 



